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Forster-Dexter theory for excitation energy transfer is generalized for the account of short time 
nonequilibrium kinetics due to the nonstationary bath relaxation. The final rate expression is 
presented as a spectral overlap between the time dependent stimulated emission and the stationary 
absorption profiles, which allows experimental determination of the time dependent rate. For a 
harmonic oscillator bath model, an explicit rate expression is derived and model calculations are 
performed in order to examine the dependence of the nonequilibrium kinetics on the excitation-bath 
coupling strength and the temperature. Relevance of the present theory with recent experimental 
findings and possible future theoretical directions are discussed. 



I. INTRODUCTION 



Excitation energy transfer (EET) is ubiquitous in 
photo sensitive materials Jp|-[l2| and is one of the key steps 
in photosynthesis j7 13-|lq|, Since the seminal works of 
Forster |],| and Dexter (FD) §, their rate expressions 
have been confirmed by numerous experiments and have 
played fundamental roles in understanding various lu- 
minescence phenomena P~p2[. In these advances, the 
spectral overlap expression of Forster jl],| , which allows 
identification of the reaction rate without any resort to 
a model Hamiltonian, has been essential. 

The rate expressions of FD are applications of the 
Fermi's golden rule (FGR), which rely on the smallness 
of the resonance interactions. The spectral overlap ex- 
pression of Forster derives from the additional simplifi- 
cation that bath modes coupled to the energy donor and 
acceptor are independent of each other. If there are com- 
mon bath modes Jl?j| , such a simple expression is not in 
general valid. However, within the harmonic oscillator 
bath model, rigorous extensions of the FD theory can 
be made [|16|-[l9| with the use of small polaron transfor- 
mation, and can be further generalized for the study of 
long-time dynamics based on the Redfield-type equation 
|ji7f and for the understanding of exciton transport in 
molecular crystals [^9|-^T|. In recent years, theoretical 
advances to cope with new experiments have been made, 
such as microscopic consideration of the medium effect 
|^2|,^3|, generalization of the FD rate expression for dis- 
ordered multichromorphic systems pl|-|30||, and unified 
theories covering up to the intermediate and strong cou- 
pling regimes [JjlJ-Q. 

Although often not noticed, the assumption of station- 
arity is implicit in the FD theory. The FGR is valid only 
when the initial density operator commutes with the ze- 
roth order Hamiltonian and the bath time scale is much 
shorter than that of the electronic transition. For EET 
processes occurring in nanosecond or longer time scale, 



one may safely assume that the bath modes have already 
relaxed and become equilibrated with the excited donor 
before the energy transfer takes place, unless there are 
ultra-slow modes or spin interactions with comparable 
time scales. Application of the FGR, thus the FD theory, 
can be justified for this case. However, EET in general is 
a nonstationary process where nonequilibrium relaxation 
of the nuclear degrees of freedom occurs during and after 
the electronic excitation. For fast EET processes occur- 
ring in a time scale comparable to the bath relaxation 
time, the reaction kinetics predicted by the FD theory 
may not be accurate enough. 

The importance of the nonequilibrium effects for fast 
EET was in fact recognized long time ago, and was named 
as hot transfer pB|-p8|. Due to experimental limita- 
tions, however, the earlier experiments and relevant theo- 
ries were concerned with the frequency domain situation 
where the reaction rate at issue is the stationary long 
time limit in the presence of an excitation field tuned for 
the hot transfer |}5 36 . Time dependent pump-probe 
situation was considered later by Sumi [^7|j3^ 1, who for- 
mulated a time dependent generalization of the FD the- 
ory based on a nonequilibrium golden rule approximation 

With the advance of ultrafast spectroscopy, it has 
become possible to induce the electronic excitation in 
the femtosecond scale. Experiments ||,|ll|,|lj,^9|-^] per- 
formed in this manner reveal systems where the EET 
rate is comparable to the vibrational relaxation rate of 
some modes. A typical example for this is the EET 
from B800 to B850 in the light harvesting complex 2 
(LH2) iQ, which is known to occur in about 1 ps. Ev- 
idences for fast EET were found in other systems also 
such as the conjugated polymer p], dendrimers B], and 
the photosynthetic reaction center ElJ. In general, the 
microscopic mechanisms of the EET in these systems are 
quite complicated, and considerations of multipolar tran- 
sitions, multichromorphic effects, and disorder may be 
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important. However, simple comparison of time scales 
indicates that examination of the nonequilibrium effects 
cannot be overlooked. The theory of Sumi |57 38 and re- 
cent theories |}2] |m| on intermediate and strong coupling 
seem suitable for these considerations. However, the na- 
ture of the approximation involved in the bath relaxation 
dynamics, upon which the nonequilibrium kinetics is sen- 
sitive, is not clear in these approaches. 

In the present work, we provide a straightforward and 
rigorous nonequilibrium generalization of the FD theory. 
The procedure is to go through the same perturbation 
theory as in deriving the FGR, but starting from the 
nonstationary initial states and considering full time de- 
pendences. Our analysis is limited to the usual pertur- 
bation regime of the FD theory. We assume that the 
resonance interaction is small enough to validate second 
order perturbation theory and that the reaction is virtu- 
ally irreversible due to either energetic or entropic reason. 
However, the excitation bath coupling is treated rigor- 
ously. The present theory is analogous to the nonequilib- 
rium generalization of the electron transfer ji^] , but an 
important distinction is that we provide a spectral over- 
lap expression valid for arbitrary bath without common 
modes, a typical situation for the EET. 

Our theory can be considered as the pump-probe ver- 
sion of the hot transfer rate theories (3^j3^]. More im- 
portantly, our spectral overlap expression brings a con- 
nection between the time dependent reaction rate and 
the modern ultrafast spectroscopy experiment, which al- 
lows direct determination of the reaction rate or enables 
experimental confirmation whether the usual assumption 
of the EET kinetics is valid. In addition, we provide cal- 
culations of the time dependent rate for the model of the 
harmonic oscillator bath, which illustrate some features 
of the nonequilibrium EET kinetics. 

The sections are organized as follows. In Sec. IIA, we 
present the main formalism and provide the spectral over- 
lap expression valid for general bath Hamiltonian without 
common modes. Section IIB provides a complementary 
result of an explicit rate expression for the harmonic os- 
cillator bath. In Sec. Ill, model calculations are made. 
Sec. IV concludes with summaries and the relevance of 
the present theory with recent experiments. 



excited to singlet states and do not consider multiexci- 
ton states. Therefore, the three electronic states \g), \D), 
and \A) form a complete set of electronic states. All the 
rest of the dynamic degrees of freedom such as molecular 
vibrations and solvation coordinates are included in the 
bath. The bath Hamiltonian is defined as that in the 
ground state \g) and is assumed to be H b — H bo + H bA , 
where the subscripts of D and A denote the components 
coupled to chromophores D and A, respectively. That is, 
in the present work, the effect of common modes is disre- 
garded. For EET processes in a medium where phonons 
and vibrations are localized, this assumption seems rea- 
sonable. 

Since we consider the excitation dynamics during a 
time much shorter than the lifetime of the excited state, 
we neglect the spontaneous decay channels of the excited 
states. Adopting a second quantization notation where 
\g) is treated as if the vacuum state, the zeroth order 
Hamiltonian describing the interaction-free dynamics can 
be written as 



(2.1) 



where e D and e A are excitation energies of D and A. 

The resonance interaction between \D) and \A) is rep- 
resented by 



H da = J{a 1 D a A + a\a D ) 



(2.2) 



where J is a function of the position vectors and the 
transition dipoles of D and A. The functional form of 
J depends on the mechanism of EET. For dipole-dipole 
interaction, it varies as the inverse third power of the 
distance between D and A. For exchange interaction, it 
is an exponential function of the distance. Here we do 
not specify the detailed mechanism, but simply assume 
that it is small enough to warrant a perturbation analysis 
and does not have any dependence on the bath operators 
(vibrational degrees of freedom), so called Condon ap- 
proximation. 

The excitation bath coupling is assumed to be as fol- 
lows: 



H eb = B D a D + B A a\ a A , 



(2.3) 



II. THEORY 

A. General bath without common modes 

The system consists of two distinctive chromophores, 
donor (D) and acceptor (A). The state where both D 
and A are in their ground electronic states is denoted as 
\g). The state where D is excited while A remains in the 
ground state is denoted as \D) = al\g), where is the 



corresponding creation operator, and \A) 



7} is de- 



fined similarly. We assume that both chromophores are 



where B D and B A are bath operators coupled to \D) and 
\A) respectively. These operators and the bath Hamilto- 
nian can be arbitrary except for the following condition: 



[Hb D ,H bA ] — [H bD1 B A ] — [H bA ,B D ] — 



(2.4) 



which implies that the bath modes coupled to \D) and 
\A) are independent of each other. This assumption can 
be justified if the chromophores D and A are far enough 
apart from each other and the major nuclear modes cou- 
pled to excitations are localized to either D or A, which 
can be consistent with the assumption of smallness of J. 
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Summing up Eqs. (2.1)-( j2.3| ), the total Hamiltonian 
describing the dynamics of the chromophores, in the sin- 
gle excitation manifold, and the bath is given by 



H = H Q + H DA + H, 



< b 



(2.5) 



For t < 0, the chromophores are assumed to be in the 
state \g) and thus the total Hamiltonian is equal to H b . 
The bath is assumed to be in the canonical equilibrium 
of H b during this period. At time zero, an impulsive 
pulse selectively excites D. Here we approximate this as 
a delta pulse in time. Then, the density operator at time 
zero, right after the irradiation, is given by 



p{Q) = \D){D\e-^/Z b 



(2.6) 



where (3 = l/k B T and Z b = Tn{e~ 0H "}. 

The probability at time t for the excitation to be found 
at A is given by 

P A (t) =Tr b UA\e- im / h p(0)e mt/h \A)\ . (2.7) 

For short enough time compared to H/J, a perturbation 
expansion of this with respect to Hda can be made. In- 
serting the first order approximation of the time evolu- 
tion operator e ~ lHt / h and the complex conjugate into Eq. 
(2.7), we obtain the following expression valid up to the 
second order of J: 



Pa® 



ft Jn 



dt' / e i{*A-*DW -*')/* 



x— Tr b 
Z b 



( e i(H b +B A )(t'-t")/h e -i(H b +B D )t'/H 



(2.8) 



The time dependent EET rate is then defined as the 
derivative of this acceptor probability as follows: 

k(t) = ±-P A {t) 



dt 
2J 2 



Re 



f dt'e l( - e °- e A )t,/h —Tr b L i{Hb+B o) t/h 
o Z b I 

x e -i(H b +B A )t'/h e -i(H b +B D )(t-t') e -pH b ^ 



where Z bA = Tr bA {e~ fm ^} and Z bo = Tr bD {e- pH "o}. 



(2.9) 



Under the assumption of Eq. (2.4), the trace over the 
bath degrees of freedom comprising Hb D and H bA can be 
decoupled from each other in the following way: 



2J 2 
k(t) = -^-Rc 

V ; h 2 

x^~Tr bA 



x^Tr b 

^br, 



^ e iH bA t'/h e -i(H„ A +B A )t'/h e - l 3H„ A | 
d ^ H b D +B D )t/h e -iH bD t' /h 



x e 



-i(H bD +B D )(t-t')/n -0H bl 



} 



(2.10) 



Equation (2.10) is the nonequilibrium generalization of 
the FD rate, but expressed in the time domain. As has 
been outlined in the introduction, this result has been ob- 
tained by applying second order perturbation theory to 
the time dependent density operator with the nonstation- 
ary initial condition of Eq. fl2.6|). The decoupled form of 



Eq. (2.10) makes it possible to express the reaction rate 
as an overlap of frequency domain spectral profiles of in- 
dependent donor and acceptor, as will be shown later. 
Before going through this procedure, it is meaningful to 
clarify the difference between the present result and the 
FD rate expression. Two additional approximations are 
necessary. 

The first is the assumption of stationarity, equivalent 
to the following replacement: 



-z(H bD +B D )(t-t')/h -0H bD i(H bD +B D )(t-t')/h 



z h 



E_ f ,-P(H bn +B D ) 



(2.11) 



where Z' b — Tr bD {e~^^ Hb D +B o)^. This approximation 
implies that the nuclear dynamics on the excited donor 
potential energy surface is ergodic in a time scale shorter 
than that of the EET transfer kinetics. We define the 
reaction rate based on this fully relaxed density operator 
as 



2J 2 

k r (t) = — ^-Rc 

y ' h 2 



dt'e i{e °- e ^ t '/ h 



hA ^ 



/h-i(H bA +B A )t'/h-0H bA 



Z b A 



—Tr bD I e i(H bD +B n )t' /h e -iH bD t' /h 

7j bD ^ 



x e 



(2.12) 



The derivation of this expression can be m ade 
by rewriting e i(H bo +B D )t/H in Eq ^ ( |2. ldj ) as 

e i(H bD + B D )(t-t')/h e i(H bD + B D )t'/h^ using tnc cyc lJ c sym _ 

metry of the trace operation for the first term, and then 
finally imposing the replacement of Eq. (2.11). 

The second is the infinite time approximation. That 
is, the FD rate corresponds to the following limit: 



k r 



k r (oo) 



(2.13) 



This approximation is valid if there is time scale separa- 
tion between the bath dynamics and the EET kinetics. 

For practical purposes of evaluating the reaction rate 
for real systems, it is important to find the explicit ex- 
pression for Eq. ( p. 10 ) in the frequency domain as an 
overlap of independent spectral profiles of the donor and 
of the acceptor. Considering the fact that the term in- 
volving the acceptor is identical in Eqs. (2.1C) and (2.12), 
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one can expect that the nonequilibrium EET rate also in- 
volves the stationary absorption profile of the acceptor. 
For this purpose, we define the absorption profile of the 
acceptor as 



I A {tu) = \fi A -e| 2 / 

J —i 



dt 



t/h_L TrbA LiH bA t/h 

ZbA L 

} , (2-14) 



reaction rate should be understood as the exponential 
decay rate of 1 — Pa(^)- The justification for this comes 
from the Redfield-type equation Q and the assumption 
that the bath relaxation has been completed. A reason- 
able way of combining these two limits is to exponentiate 
the time integration of k(t) as follows: 



x e -i(H bA + B A )t/Ti e -f3H hA 



where fx a is the transition dipole of the acceptor and e 
is the polarization vector of the radiation. Inserting the 
inverse transform of this into Eq. ( [2.10 ) , 



J 2 

k(t) = — 5, TT 2 



dto Ia(w) 



x Re 



,)t/n 



[ dt'e- luJt ' +t,i n t ' /h —Tr bD [ e l(HbD+B ° ] 

JO ZbD L 

x e -iH bn t'/h e -i(H b n+B D )(t-t')/h e -0H b n\~ _ (2.15) 

The time dependent part involving the donor can be ex- 
pressed as the stimulated emission profile in a pump- 
probe experiment. In Appendix A, we derive a time de- 
pendent stimulated emission profile of the donor which is 
subject to a stationary field after being excited by a delta 
pulse. Inserting Eq. (A5) into Eq. (2.15), the frequency 



domain expression of the EET rate is given by 



k(t) 



J 2 



2irh \fiA ■ e\ 2 \nn • e 



ai2 



dujlA(u)E[)(t,uj) . 

(2.16) 



In the limit where t — > oo, this expression becomes equiv- 
alent to the Forster's spectral overlap expression as long 
as Ed(oo,u>) is equal to the spontaneous emission pro- 
file of the excited donor except for a normalization factor 
and the univ ersal frequency dependent scaling function. 

Equation ( [2.16| ) is the central result of the present pa- 
per. It is the pump-probe version of the hot transfer 
rate |35],[36| and generalizes the spectral overlap expres- 
sion of Forster for fast EET processes. With the modern 
development of ultrafast spectroscopy, determination of 
k(t) and Ejj (t, ui) can be made even in femtosecond scale. 
With the advances in experimental techniques of altering 
chromophores by chemical or biological manipulations, 
independent determinations of I A (lo) and ED{t,ui) can 
be done at the same condition as that of kit) for a broad 
range of systems. If these practical issues are settled and 
if the system satisfies the requirements for applying the 
perturbation theory, Eq. (2.16) should hold as long as 
the effect of common modes is insignificant. 

In Eq. (2.S), we have defined the reaction rate as the 
time derivative of P A (t)- Due to the use of perturba- 
tion theory, such a definition gives a valid result only for 
J Q dt'kit') « 1. In the longer time limit when the pop- 
ulation transfer has occurred significantly, instead, the 



P A (t) w 1 - exp 



dt' k(t') 



(2.17) 



where it has been assumed that the transfer from D to A 
is irreversible due to either energetic or entropic reason. 
Equation (2.17) is not based on a rigorous derivation, 
and conditions when such approximation is valid need to 
be clarified based on a more rigorous approach, which 
will be done elsewhere. However, for the purpose of un- 
derstanding the qualitative aspect of EET kinetics in the 
nonequilibrium situation, which is the main purpose of 



the present paper, the expression of Eq. (2.17) is useful 



B. Linearly coupled harmonic oscillator bath 

For the simple case where the bath consists of inde- 
pendent harmonic oscillators, an explicit expression can 
be found for k(t). Assume that the bath Hamiltonian is 
given by 



(2.18) 



and the chromophore-bath interaction is linear in the 
bath coordinate as follows: 

H e b = ^huj n (b n + bl l )(g nD a^ D a D + g nA a\a A ) , (2.19) 

n 



where g n D9nA — and thus the condition of Eq. (2.4) 
is satisfied. For an explicit calculation, it is convenient 



to start from Eq. (2.10). The integrand of the reaction 



rate involves trace of the product of the propagators for 
displaced harmonic oscillators. Each trace over the ac- 
ceptor bath and the donor bath can be done explicitly, 
and the resulting expression for the reaction rate can be 
written as 



2 J 2 
kit) = -^-Rc 

h 2 







X e 2 *^-^ al D ( s ' m (^rzt)-sin(uj n (t-t')) 

X e~ 2 T,Jaln+9l A ) coth(/3^„/2) sm 2 (u n t'/2) 



(2.20) 



where 



e D(A) = e D(A) 



, hUn9n D {A) 



£d(A) 



- A 



D(A) 



(2.21) 
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where A u and are reorganization energies of the donor 
and acceptor baths. For the harmonic oscillator bath 
model, in fact, the reaction rate can be calculated ex- 
plicitly even for more general situation where there are 
common modes coupling both the donor and acceptor. 
In Appendix B, we derive the general expression employ- 
ing the small polaron transformation. Equation ( 2.20| ) 
can alternatively be obtained from Eq. (BIO) by letting 
g n Dg n A = 0. 

Define the following spectral densities: 



n 



(2.22) 
(2.23) 



Inserting these definitions into Eq. (2.20), the nonequi- 
librium EET rate can be expressed as 



k(ty- 



2J 2 



-Re 



xe 



xe 



-i duj(r}D (u)+r]A (^-0) sin(ujt' ) / uj 2 
2i dwqjj (oj) (sin(aji) — sin(u;(i — t'))) /u 2 

> duj(rin(u)+r)A(u)) coth(0Hu> /2) sin 2 (wi'/2)/w 2 



(2.24) 

This is the main result of the present subsection. The 
expre ssion for k r (t), which can be calculated from Eq. 
( 2.12 ) through a similar procedure is given by 



2 J 2 

k r (t) = — o" Re 

h 



dt' e i(l n- l A )t,/h 



(I 



^_ — i f a duj(r]r> (uj)-\-r]A (<*>)) sin(ujt' ) / u> 2 

-2 f 00 du>(r)n(u)+VA(u)) coth(/3w/2) sin 2 {tot' /2) /u 



(2.25) 

Due to the nonstationary initial condition, the integrand 
in Eq. ( 2.24 ) has an additional dependence on t, which 
becomes clear when we compare the expression with that 
of k r (t). In general, further simplifications of the rate 
expressions given by Eqs. ( 2.24 ) and ( 2.25 ) cannot be 
made. However, this does not pose a practical difficulty 
in calculating the reaction rate because direct numerical 
integrations of these expressions can be done quite easily. 
In the next section, we implement these calculations for 
a simple model spectral density. 

Before concluding this section, we examine one impor- 
tant limit where the stationary phase approximation is 
possible. In the strong excitation-bath coupling or high 
temperature limit, the dominant contribution of the in- 
tegration comes from small t' region. Expanding all the 
functions of t' in the exponent up to the second order, 



k s (t) 



2J< 



Re 



2J^Re 

x f dt'e 
Jo 



fai e i(6e-X T +2C(t))t' /h-(D(fj)/2-iS(t))t' 2 



-(8e-\ T +2C(t)) 2 /(2h 2 (D(f3)-2iS(t))) 



, -( D(0)/2-iSV)) (t'-i ^fXfsUl ) 



where 

Jo ^ 
S(t) = / cLor]£)(uj) sm(uit) , 



A_d + A^ 



C(t) = h 



au> cos{Lut) 



dui(r]D{u) + t]a{u)) coth( 



(2.26) 

(2.27) 
(2.28) 
(2.29) 
(2.30) 



Although Eq. (2.26) involves a simpler integration than 
Eq. (2.24), the complex Gaussian integration for finite 
t does not convey a clear physical picture. If \S{t)\ « 
D(f3) and in the long enough time limit, one can disregard 
the imaginary term S(t) and approximate the integration 
in Eq. (2.26) as that from to oo. The rate expression 
under this situation can be simplified as 



k S g(t) 



J 



I 2tt 



^-(Se-\ T +2C(t)) 2 /2h 2 D{0) 



(2.31) 



Eq. (2.24) can be approximated as 



This approximation is valid only for long time regime 
and for strong enough excitation-bath coupling. The 
maximum rate based on this expression is obtained for 
6e = At — 2C(t). Since C(oo) = for most realistic 
spectral densities (cf. sub-Ohmic, C(oo) ^ 0), this re- 
lation implies that the optimum energy difference in the 
long time limit is given by Se = At- However, the above 
expression suggests that the reaction rate for a nonop- 
timum 8e can also be substantial during the transient 
period when C{t) changes. If J is not small compared 
to the bath relaxation rate, the contribution of this tran- 
sient term cannot be neglected. A similar observation 
was made in the theory of noncquilibrium electron trans- 
fer reaction IM. 



III. MODEL CALCULATIONS 

Assume that the bath modes coupled to D and A have 
the same spectral profile, while their magnitudes of cou- 
pling can be different. We consider a model where the 
number of modes increases linearly for small lu and has 
an exponential tail for large uj. The corresponding spec- 
tral densities, defined by Eqs. (2.22) and (2.23), thus 
have the following form: 
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(3.1) 



where a D and a A determine the coupling strengths of 
\D) and \A) to their respective bath modes. u> c is the 
cutoff frequency which determines the spectral range of 
the bath. Due to the common value of this cutoff fre- 
quency, the bath relaxation dynamics has a single time 
scale and the reaction rate exhibits a simple scaling be- 
havior. In Appendix C, we provide a detailed form of the 
reaction rate in the units where K = 1. As can be seen 
from Eq. (CI), all the EET kinetics for different ui c can 



be deduced from the following dimcnsionless rate: 



T = LO r t 



(3.2) 



For the case where the reaction rate is time dependent, 
the instantaneous value of the reaction rate cannot be a 
clear indication of the efficiency of the energy transfer. 
For this reason, we also provide the following dimension- 
less quantity: 



p( T ) = / <It'k{t') 
Jo 



(3.3) 



According to Eq. fl2.17| ), P A (t) « l-exp{-( J/uj c ) 2 p(t)}, 
the population of \A) at time t — t/uj c . 

For the reaction rate k r (t) given by Eq. (2.25), sim- 



ilar quantities can be defined and these are denoted as 
Kf(r) and p r {r). The expression for K r (r) can be ob- 
tained from that of k(t) simply disregarding the third 
term in Eq. (|C2| ) . Equations (C5)-(|C7|) in combination 
with Eqs. ( 2.26| ) and ( [2.31 ) provide explicit expressions 
for k s (t) and k sg (t) . The dimcnsionless versions of these 
rates are denoted as k s (t) and K sg (r) , and their cumula- 
tive integrations are represented by p s (t) and p sff (r). 



A. Zero temperature limit 

We first performed calculations for = a a = 10, 
a strong coupling limit. The values of the renormalized 
energy difference Se were chosen to be At — Xd , A y, an d 
At + A_d , which were motivated by the form of Eq. (2.31). 



These three cases are analogues of the normal, activation- 
less, and inverted regimes in the electron transfer kinet- 
ics. Figure 1 shows the rates. The nonequilibrium effect 
is shown to be important until about r « 2. This is in 
contrast to the behavior of K r (r), which approaches the 
long time limit in about r « 0.2. During the nonsta- 
tionary period, k(t) goes through a maximum before ap- 
proaching the long time limit in a smooth fashion. The 
maximum value appears earlier for smaller value of <5e, 
for which it takes shorter for the bath to relax into the 
resonance condition. In the classical limit, this can be un- 
derstood in terms of the average time for the wavepacket 
in the excited state to go through the crossing point. 



The two approximations of k s (t) and k S9 (t) repro- 
duce the long time limit quite well for the present case of 
strong coupling limit. In the nonstationary regime, they 
exhibit a qualitative behavior similar to that of k(t), but 
the quantitative agreement is not satisfactory in the re- 
gion of 1 < r < 2. This indicates that, in this interme- 
diate time regime, the coherence of the bath still plays a 
substantial role and cannot be well approximated by the 
stationary phase integration. 

In Fig. 2, p(r) given by Eq. ( |Q| ) and analogous 
quantities for other approximations are provided. For 
Se = Xt — Xd, the nonequilibrium contribution results 
in larger population transfer during the short initial pe- 
riod, but not in the longer time period. For Se = Xt, the 
nonequilibrium effect always leads to smaller population 
transfer. For Se = Xt + Xd, the nonequilibrium popula- 
tion transfer is smaller than the fully relaxed one until 
about t ~ 1, but it becomes larger during the longer time 
as the delayed response of the bath brings a resonance 
condition. The different kinetics during the nonstation- 
ary period is seen to result in finite differences in the net 
population transfer in the long time regime. The values 
of Ps(t) and p sg {r) are seen to overestimate p(r) when 
Se = Xt — Xd and underestimates it for Se = Xt + Ad- 
The stationary phase approximation seems to work rela- 
tively well for the case of Se = Xt- 

Some of the patterns observed from the above calcula- 
tions may be specific for the model of the bath and the 
temperature. However, the features that the nonequilib- 
rium effect is substantial until about up to r w 2 and 
that the maximum peak of the reaction rate appears ear- 
lier for smaller Se are expected to hold quite generally. 
Due to the assumption of strong coupling, the approx- 
imate expressions of k s (t) and K sg (r) agree quite well 
with k(t), although there is still a difference in the net 
population. Such agreement cannot be seen for the weak 
coupling case considered next, an expected result consid- 
ering the assumptions involved in k s (t) and K sg (r). 

Calculations were performed for ac = 1 and a a = 1, 
a weak coupling limit. Figure 3 shows the dimcnsionless 
reaction rates and Fig. 4 their cumulative integrations. 
Unlike the case of strong coupling, the difference between 
k(t) and n r (r) is not appreciable except for some phase 
difference in their oscillations. The resulting time inte- 
grations, p(t) and p t {t), are quite close as can be seen 
from Fig. 4. Especially for the two cases of Se = Xt 
and Se = Xt + Xd, the two curves almost overlap with 
each other. The approximations of k s (t) and k s3 {t) are 
not expected to work well for the present weak coupling 
regime, but it is meaningful to examine their qualitative 
nature. In Fig. 3 , it is shown that k s (t) and k S9 (t) do 
not reproduce the oscillatory behavior and deviate from 
k(t) systematically, except for the case of Se = Xt- The 
resulting values of p s (r) and p sg {T) underestimate p(r) 
for Se = Xt — Xd and overestimate it for Se = Xt + Xd- 
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These are opposite to the trends observed for strong cou- 
pling limit. In conclusion, for the present case of weak 
coupling, the nonequilibrium effect is unimportant, and 
both k(t) and K r (r) rise to their steady state limit in a 
similar fashion. The oscillations persistent in both k(t) 
and k,.(t) in Fig. 3 originate from e j( - £ o~ e ^- )t in Eqs. 
( [2.24 ) and ( 2.25| ), which remains significant even in the 
long time limit due to the weak excitation-bath coupling. 
However, even though all the assumptions of the present 
section hold, actual observation of these oscillations in 
a real system is unlikely due to the energy uncertainty 
related dephasing. 



B. Temperature effect 

As the temperature increases, the phonon side band 
becomes broader in the spectral overlap expression of 
Eq. (2.16). The relaxation and dephasing of the bath 
become faster also. It is interesting to examine how these 
changes affect the nonequilibrium kinetics. For this pur- 
pose, we compare only k(t) and K r (r). The temperature 
dependence was accounted for by using an approximate 
expression for 7(r), in Eq. (|Cl|), as follows: 



7(t) » (a D + a A ) 
+ 



2r 2 



1 + r 2 (1 + T 2 ) 2 

16r 2 



/3w c (/3w c + 2)((/?w c + 2) 2 + 4r 2 ) 



(3.4) 



This is based on the approximation of coth(a;/2) 
1 + (2/x)e~ x / 2 , which reproduces the proper low and 
high temperature limits. Although the approximation 
is rather crude in the region of (3uj c ~ 1, the overall per- 
formance is good enough to produce a correct qualitative 
trend. 

In Fig. 5, we provide results for an — ola = 10, 
the strong coupling limit, and for Se — Xt — Xd- The 
zero temperature result is shown in the top panel of Fig. 
1. As the temperature increases, the nonequilibrium ef- 
fects diminish. However, even for very high tempera- 
ture limit of j3us c = 0.1, there is a noticeable difference 
between k(t) and K r (r). Only in the highest tempera- 
ture limit of [3uj c = 0.01, k(t) and K r (r) show an agree- 
ment. This indicates that for strongly coupled systems, 
the nonequilibrium effects persist even up to very high 
temperature. Next, we considered the weak coupling case 
of old = a a = 1 with the same choice of Se — At — Xp. 
The results are shown in Fig. 6. In the zero temperature 
limit, the nonequilibrium effect was not so significant for 
this weak coupling case. This fact remains the same even 
for the higher temperature. 



IV. DISCUSSION 

The main focus of the present paper was the nonequi- 
librium bath relaxation effect. For this reason, we con- 
sidered only the simplest case where donor and acceptor 
each have one excited state. However, as recent stud- 
ies on various systems indicate |7|,[l3|-|l5|] , the existence of 
multiple chromophores is a general situation rather than 
an exception. Therefore, one in general must consider 
the case where the donor and the acceptor respectively 
consist of multiple exciton states. Given that the inter- 
excitonic relaxation can be disregarded or approximated 
in a self consistent way, the generalization of the present 
theory can be done straightforwardly as are those of the 
FD theory p4]-|30|l . In many cases, this type of simple ex- 
tension may indeed capture the main aspect of the EET 
between multichromorphic donor and acceptor. However, 
if some of the excitonic levels are closely spaced or if there 
is degeneracy due to symmetry, more careful theoretical 
analysis is necessary. 

In a heterogeneous environment or a complex system, 
disorder plays an important role and its explicit consid- 
eration becomes necessary for a proper understanding of 
the system. Due to the averaging over the static disor- 
der involved in this case, our spectral overlap expression 
cannot be used directly for the observables of an ensem- 
ble experiment. However, given that the distribution of 
the static disorder is well characterized by other means, 
Monte Carlo simulation can provide an indirect way of 
determining the nonequilibrium population transfer. The 
study of the interplay between the heterogeneity and the 
nonequilibrium effect in this way can bring new insights 
into the fast EET kinetics of complex systems. 

In the LH2 of purple bacteria, the EET from B800 to 
B850 occurs in about lps and the rate is not so sensi- 
tive to temperature fl4|| . This rate is much faster than 
that obtained by the FD theory, but it has recently been 
shown that consideration of the multichromorphic nature 
of the B850 and the disorder account for much of the dis- 
crepancy pjj-f??! . However, the theoretical value is still 
somewhat smaller than the experimental one [^7) . Many 
explanations are possible for this, and the nonequilibrium 
effect is one such possibility. 

According to our model calculations, the nonequilib- 
rium short time kinetics shows complicated time depen- 
dent behavior and the reaction rate in this regime is less 
dependent on the spectral overlap between the stationary 
emission and absorption profiles. Some of these features 
can be seen in recent sub-picosecond pump-probe exper- 
iments ||,[ll],|4]]. For example, the reaction rate is shown 
to be relatively insensitive to the change of the spectral 
profile for the EET in the photosynthetic reaction center 

& 

Before applying the present theory to a given system, 
it is always important to examine whether the nature 
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of the system is consistent with the assumption of irre- 
versibility and the use of second order approximation. It 
is expected that our results can be applied as long as 5e is 
sufficiently large and the bath relaxation of the acceptor 
is fast enough. However, in order to address these issues 
more systematically, it is necessary to formulate a the- 
ory that can account for electronic coherence. For this 
purpose, one may adopt the formalism of the generalized 
master equation E3| or consider in the framework of the 
Rccffield-type equation p7[. 
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APPENDIX A: TIME DEPENDENT 
STIMULATED EMISSION PROFILE OF THE 
DONOR 

Consider a system consisting of the donor and its own 
bath. The relevant Hamiltonian in the single excitation 
manifold is given by 



HD=e D a ] D a D +B D a ] D a D +H bD . 



(Al) 



Assume that the donor is excited at time zero by a delta 
pulse. Then, the initial density operator at t = is given 
by 



p (0) = \D)(D\^e-^o 



(A2) 



Right after the pulse excitation, a stationary field with 
a fixed frequency lu is turned on. The time dependent 
Hamiltonian governing the dynamics for t > 0, assuming 
unit field strength, is given by 

H(t) = H D + \^ D -e\(e-^ t \D)(g\ + e iut \g)(D\) , (A3) 

where rotating wave approximation was used. The prob- 
ability for the donor to emit a photon and go to the 
ground state, induced by the stationary field, is given by 



P g (t) = \fx D ■ e\ 2 dt' dt" e Mt'-t")-ie D (t'-t")/n 
Jo Jq 

x _LjY, / e -iH bD (t-t')/h -i(B D +H„ D )t' /K 

z hn bD \ 



J bD 
x e 



E D (t,tu) 



d 
It 



2\n D -e| 2 Re 



dt 'e~ iuJt ' +i£ c t ' l h 



:-^—Tr bo I e i(Hb D +Bn)t/H e -iH bn t ' /h 

Z b D ° L 

-i(H bD +B D )(t-t')/h e -/3H bD Y _ ( A5 ) 



x e 



APPENDIX B: RATE EXPRESSION FOR A 
GENERAL HARMONIC OSCILLATOR BATH 

For a general harmonic oscillator bath model, we de- 
rive the expression of the nonequilibrium reaction rate 
given by Eq. based on the small polaron transfor- 

mation. For this purpose, first we introduce the following 
generator of small polaron transformation: 

S = - ^2{bn - b\ l ){g nD a) D a D + g nA a\a A ) . (Bl) 

n 

Then, it is straightforward to show that 



H = e S He S = e D a^ D a D + l A a\a A + H b + H 



DA 



where 



H DA = Je s {a [ D a A + a\a D )e~ s 

= J{ey A (J D a A +e\6 D a\a D ) 



where 



'D(A) 



(B2) 



(B3) 



(B4) 



Inserting 1 = e s e s between every two operators in Eq. 
.91), one can show that 



2 J 2 

k(t) = -^Rc 

h 2 



x—Tr b {e lHbt 'l n &e- lHbt ' /h Qp d {t - t')\ 



(B5) 



where 0t = 9l 



with 



d w a, 

p d (t) = e- lHbt ' h e^ D e- f3Hb e D e lHbt / h 



b n D(t) =b n — g n De~ 



(B6) 
(B7) 



l3H„ De i(B D +H bD )t" /h e iH bD (t-t")/h^ ^ (- A4 j 



where weak field approximation was used. The time de- 
pendent stimulated emission profile is defined as the time 
derivative of this probability as follows: 



In Eq. (|B5|), 

e iH b t/hQ^ e -iH b t/h 

= exp < -^.g„A0nf 

I n 



(B8) 
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where g n A = g n D — 9nA- Inserting Eqs. (B6)-(B8) into 
Eq. ©, 



k(t) 



2J 2 



Re 



[ dt'e^n-^'^e-^n'- 
Jo 



In the main text, we consider the zero temperature limit 
first and study the finite temperature effect using a sim- 
ple interpolation formula. Similarly, the quantities enter- 



i(u„t') 



x — Tr b { e 



x e 



where the time dependences of b n D{t — t') and b' nD (t — t') 
have been omitted, but this does not make any difference 
in the result because of the trace operation. Disentan- 
gling the summation in the exponent and the n ev aluating 
the trace over the bath, one can reduce Eq. (B9) into the 
following integral expression: 



k(t) 



2J 2 



-Re 



" 2 2 coth(/8fiui n /2) sin 2 (w 7l t72)" 



xe 



(BIO) 



For the case where there is no common mode, 
9nA9nD = gin and gl A = g 2 nD + g 2 nA and the expres- 
sion of Eq. (2.20) is reproduced. 



APPENDIX C: RATE EXPRESSIONS FOR THE 
MODEL SPECTRAL DENSITIES OF SEC. Ill 



For the spectral densities given by Eq. (3.1), the in- 
tegrations within the exponents of Eq. ( [2.24 ) can be 
performed explicitly Adopting the units where h = 1, 



the reaction rate of Eq. (2.24) can be expressed as 



2 J 2 

hit) = —Re 



fal e i0(r,r')-7(T"') 



(CI) 



where r = uj c t, t' = uj c t' , and 



- 2(a D - at; 



2{a D +a A y 
(1 



r /2\2 



-4a r 



(l + (r-r') 2 ) 2 (l + T 2 f 



(C2) 



In Eq. (CI), the function j(t') comes from the integra- 
tion involving coth(/3u;?i/2) in Eq. ( ^.24 ). Simple ap- 
proximation for this function is possible in the limits of 
low and high temperature. Using these approximations, 



7(r)= 



(a D + a A ) [j^pi + (!+*»)* ) ' » 1 

2(a D +a A ) r 2 

— g jjppz , fJUJc « 1 • 



-l)6„ D -(e iw » t -l)6t D ) 



(B9) 



ing the approximate rate expressions of Eqs. (2.26) and 
( [2.31 ) can be calculated. These are 



X T = 2(a D + a A )ui, 
S(t) = 24a D u 2 c ^~ 



(l+T 2 ) 4 



C(t) = 2a D u 
D(J3)-- 



1 



(l+T 2 ) 



2\3 



LO r t 



T = U) r t 



6(a D + a A )uj 2 , (3l} c » 1 



4(a D +a A )u c /P , I3u c « 1 



(C4) 
(C5) 

(C6) 

(C7) 



(C3) 
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FIG. 1. Dimensionless nonequilibrium reaction rate k(t), Eq. (3.2), and analogous quantities for k r (t), k s (t), and k sg (t). 
a D = a A — 10 and temperature is zero. The top panel corresponds to Se — At — Ad, the middle panel to Se = At, and the 
bottom panel to Si = At + Ad. 
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FIG. 3. Dimensionless nonequilibrium reaction rate k(t), Eq. (3.2), and analogous quantities for k r (t), k s (t), and k ag (t). 
ol d = <x A = 1 and temperature is zero. The top panel corresponds to Se = At — Ad, the middle panel to 8e = At, and the 
bottom panel to 8e = At + Ad. 
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